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ABSTRACT 

Almost none of the r-modes ordinarily found in rotating stars exist, if the star 
and its perturbations obey the same one-parameter equation of state; and rotating 
relativistic stars with one-parameter equations of state have no pure r-modes at all, no 
modes whose limit, for a star with zero angular velocity, is a perturbation with axial 
, parity. Similarly (as we show here) rotating stars of this kind have no pure g-modes, 

On I no modes whose spherical limit is a perturbation with polar parity and vanishing 

Q I perturbed pressure and density. Where have these modes gone? 

CN ' In spherical stars of this kind, r-modes and g-modes form a degenerate zero- 

QQ ■ frequency subspace. We find that rotation splits the degeneracy to zeroth order in the 

■ star's angular velocity $1, and the resulting modes are generically hybrids, whose limit 

O . as — > is a stationary current with axial and polar parts. Lindblom and Ipser have 

I . recently found these hybrid modes in an analytic study of the Maclaurin spheroids. 

Since the hybrid modes have a rotational restoring force, they call them "rotation 
modes" or "generalized r-modes". 

Because each mode has definite parity, its axial and polar parts have alternating 
;h ' values of I. We show that each mode belongs to one of two classes, axial-led or 

polar-led, depending on whether the spherical harmonic with lowest value of / that 
contributes to its velocity field is axial or polar. We numerically compute these modes 
for slowly rotating polytropes and for Maclaurin spheroids, using a straightforward 
method that appears to be novel and robust. Timescales for the gravitational-wave 
driven instability and for viscous damping are computed using assumptions appropriate 
to neutron stars. The instability to nonaxisymmetric modes is, as expected, dominated 
by the / = m r-modes with simplest radial dependence, the only modes which retain 
their axial character in isentropic models; for relativistic isentropic stars, these / = m 
modes must also be replaced by hybrids of the kind considered here. 
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1. Introduction 

A recently discovered instability of r-modes of rotating stars (first found in numerical study 
by Andersson (1998), and analytically verified by Friedman and Morsink (1998)) has gained the 
attention of a number of authors ( Kojima 1998| ; Lindblom, Owen, and Morsink 1998; Owen et 



al. 1998; Andersson, Kokkotas and Schutz 1998; Kokkotas and Stergioulas 199^ ; Andersson 



Kokkotas and Stergioulas 1998; Lindblom and Ipser 1998; Madsen 1998; [Spruit 1999| ; Beyer and 



Kokkotas 1999 ; Lindblom et al. 1999 ). These modes have axial parity (see below) and their 
frequency is proportional to the star's angular velocity. Neutron stars that are rapidly rotating at 
birth are likely to be unstable to nonaxisymmetric perturbations driven by gravitational waves; 
estimates of growth times and viscous damping times (Lindblom et al 1998, Pwen et al. 199^ , 
Andersson et al 1998, Kokkotas and Stergioulas 1998| , [Lindblom et al. 1999| ) suggest that r-modes 



dominate the spin-down of such stars for several months, until a superfluid transition shuts off the 
instability. Unstable r-modes may thus set the upper limit on the spin of young neutron stars, and 
gravitational waves emitted during the initial spin-down might be detectable. The recent discovery 
by Marshall et al (1998) of a pulsar in the supernova remnant N157B implies the existence of a 
class of neutron stars that are rapidly rotating at birth and whose spin is plausibly limited by the 
gravitational-wave driven instability. 

Perturbations of a spherical star can be divided into two classes, axial and polar, depending 
on their behavior under parity. Where polar tensor fields on a 2-sphere can be constructed from 
the scalars Yj^ and their gradients VYJ"* (and the metric on a 2-sphere), axial fields involve the 
pseudo- vector f x VY^"^, and their behavior under parity is opposite to that of Y[^. That is, axial 
perturbations of odd I are invariant under parity, and axial perturbations with even / change sign. 
If a mode varies continuously along a sequence of equilibrium configurations that starts with a 
spherical star and continues along a path of increasing rotation, the mode will be called axial if it 
is axial for the spherical star. Its parity cannot change along the sequence, but / is well-defined 
only for modes of the spherical configuration. 

It is useful to further divide stellar perturbations into subclasses according to the physics 
dominating their behaviour. In perfect fluid stellar models, the polar parity perturbations consist 
of the f-, p- and g-modes; the f- and p-modes having pressure as their dominant restoring force 
and the g-modes having gravity as their dominant restoring force (Cowling 1941). The axial parity 
perturbations are dominated by the Coriolis force in rotating stars and were called "r-modes" 
by Papaloizou and Pringle (1978) because of their similarity to the Rossby waves of terrestrial 
meteorology. For slowly rotating stars, the r-modes have frequencies which scale linearly with 
the star's angular velocity and in the spherical limit they become time-independent convective 
currents. 

Despite the sudden interest in these modes, however, they are not yet well-understood for 
stellar models in which both the star and its perturbations are governed by a one-parameter 
equation of state, p = p{p); we shall call such stellar models isentropic, because isentropic models 
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and their adiabatic perturbations obey the same one-parameter equation of state. For stars with 
more general equations of state, the r-modes appear to be complete for perturbations that have 
axial parity. However, this is not the case for isentropic models. One finds that the only purely 
axial modes allowed in isentropic stars are the physically interesting / = m r-modes with simplest 
radial behavior ( [Papalouizou and Pring le 19781 ; Provost et al. 1978Q; | aio 1982] ; pmeyers and] 



Martens 1983 ). The disappearance of the purely axial modes with / > m occurs for the following 
reason. In spherical isentropic stars the gravitational restoring forces that give rise to the g-modes 
vanish and they, too, become time-independent convective currents with vanishing perturbed 
pressure and density. Thus, the space of zero frequency modes, which generally consists only 
of the axial r-modes, becomes larger for spherical isentropic stars to include the polar g-modes. 
This large degenerate subspace of zero-frequency modes is split by rotation to zeroth order in 
the angular velocity, and the corresponding modes of rotating isentropic stars are hybrids whose 
spherical limits are mixtures of axial and polar perturbations. These hybrid modes have already 
been found analytically for the uniform-density Maclaurin spheroids by Lindblom and Ipser (1998) 
who point out that since their dominant restoring force is the Coriolis force, it is natural to refer 
to them as rotation modes, or generalized r-modes. 

Although isentropic Newtonian stars do retain a vestigial set of purely axial modes (those 
having / = m), it appears that rotating relativistic stars of this type have no pure r-modes, no 
modes whose limit for a spherical star is purely axial (Andersson, Lockitch and Friedman, 1999). 
For nonisentropic relativistic stars, Kojima (1998) has derived an equation governing purely axial 
perturbations to lowest order in the star's angular velocity.^ For the isentropic case, however, 
we find a second, independent equation for these perturbations that appears to give inconsistent 
radial behaviour for purely axial modes (Andersson, Lockitch and Friedman, 1999). This second 
equation is obtained from an angular component of the curl of the relativistic Euler equation and 
is simply the relativistic generalization of the equations presented in Sect. Ill of this paper (Eq. 
(p^), for example). Kojima's equation is one from which polar-parity perturbations have been 
excluded, and, as the present paper makes clear, one cannot assume that axial and polar parity 
modes decouple for slowly rotating isentropic stars. Instead, we expect the Newtonian I = m 
r-modes to become discrete axial-led hybrids of the corresponding relativistic models. 

In this paper we examine the hybrid rotational modes of rotating isentropic Newtonian stars. 
We distinguish two types of modes, axial-led and polar-led, and show that every mode belongs to 
one of the two classes. We then turn to the computation of eigenfunctions and eigenfrequencies 
for modes in each class, adopting what appears to be a method that is both novel and robust. For 



^An appendix in this paper incorrectly claims that no / = m r-modes exist, based on an incorrect assumption 
about their radial behavior 

^Based on this equation, Kojima has argued that the spectrum is continuous, and his argument has been made 
precise in a recent paper of Beyer and Kokkotas (1999) (See also Kojima and Hosonuma 1999). Beyer and Kokkotas, 
however, also point out that the continuous spectrum they find may be an artifact of the fact that the imaginary 
part of the frequency vanishes in the slow-rotation limit. 
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the uniform-density Maclaurin spheroids, these modes have been found analytically by Lindblom 
and Ipser in a complementary presentation that makes certain features transparent but masks 
properties that are our primary concern. We examine the eigenfrequencies and corresponding 
eigenfunctions to lowest nontrivial order in the angular velocity $7. We then examine the 
frequencies and modes of n = 1 polytropes, finding that the structure of the modes and their 
frequencies are very similar for the polytropes and the uniform-density configurations. The 
numerical analysis is complicated by a curious linear dependence in the Euler equations, detailed 
in Appendix B. The linear dependence appears in a power series expansion of the equations about 
the origin. It may be related to difficulty other groups have encountered in searching for these 
modes. 

Finally, we examine unstable modes, computing their growth time and expected viscous 
damping time. The pure I = m = 2 r-mode retains its dominant role, but the 3 < i = m ^ 10 
r-modes and some of the fastest growing hybrids may contribute to the gravitational radiation and 
spin-down. 



2. Spherical Stars 

We consider a static spherically symmetric, self-gravitating perfect fluid described by a 
gravitational potential density p and pressure p. These satisfy an equation of state of the form 

P = Pip), (1) 
as well as the Newtonian equilibrium equations 

V„(/i + $)=0 (2) 

= AnGp, (3) 
where h is the specific enthalpy in a comoving frame. 




We are interested in the space of zero-frequency modes, the linearized time-independent 
perturbations of this static equilibrium. This zero-frequency subspace is spanned by two types of 
perturbations: (i) perturbations with Sv°- ^ and Sp = Sp = = 0, and (ii) perturbations with 
Sp, 6p and (5$ nonzero and Sv"" = 0. If one assumes that no solution to the linearized equations 
governing a static equilibrium is spurious, that each corresponds to a family of exact solutions, 
then the only solutions (ii) are spherically symmetric, joining neighboring equilibria. 

The decomposition into classes (i) and (ii) can be seen as follows. The set of equations 
satisfied by {6p,5^,5v"') are the perturbed mass conservation equation. 



6[dtp + Va{pv'')]=0, 



(5) 
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the perturbed Euler equation, 

\dt + £,)Va+ya{h-\v^ + ^)\=Q, (6) 

and the perturbed Poisson equation, J[Eq. @]. 

For a time-independent perturbation these equations take the form 

VaipSv'') = 0, (7) 
Va{5h + 5^) = 0, (8) 

and 

V2<^$ = 47rG(5/), (9) 

where 

<^/^ = ^ = ^^. (10) 
P dp p 

Because Eq. (|7|) for 5v^ decouples from Eqs. (P), (P) and (|To|) for 5$), any solution to 
Eqs. (0)-(lO) is a superposition of a solution (0,0, and a solution (5/?, 5$,0). This is the 
claimed decomposition. 

The theorem that any static self-gravitating perfect fluid is spherical implies that the solution 
((5p, 5^>,0) is spherically symmetric, to within the assumptions that the static perturbation 
equations have no spurious solutions ("linearization stability" )0. 

Thus, under the assumption of linearization stability we have shown that all stationary 
non-radial perturbations of a spherical, isentropic star have 5p = 6p = 5^ = Q and a velocity field 
5v°- that satisfies Eq. (^). 

A perturbation with axial parity has the form ( Friedman and Morsink 1998| ), 

Sv"" = Uir)e'''"'VbYrVcr, (11) 

and automatically satisfies Eq. 

A perturbation with polar parity perturbation has the form, 

^^a ^ 1^:1!:) ymya^ ^ V{r)V''Yf^; (12) 

r 

and Eq. (|^ gives a relation between W and V, 

^(rpW)-l(l + l)pV = 0. (13) 
dr 



^We are aware of a proof of this linearization stability for relativistic stars under assumptions on the equation of 
state that would not allow polytropes (Kiinzle and Savage 198C). 
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These perturbations must satisfy the boundary conditions of regularity at the center, r = 
and surface, r = R, the star. Also, the Lagrangian change in the pressure (defined in the 
next section) must vanish at the surface of the star. These boundary conditions result in the 
requirement that 

W{0) = W{R) = 0; (14) 
however, apart from this restriction, the radial functions U{r) and W{r) are undetermined. 

Thus, a spherical, isentropic, Newtonian star admits a class of zero frequency convective fluid 



motions of the forms ( pT| ) and (12). Because they are stationary, these modes do not couple to 
gravitational radiation. ^ 



3. Rotating Isentropic Stars 

We consider perturbations of an isentropic Newtonian star, rotating with uniform angular 
velocity 0,. No assumption of slow rotation will be made until we turn to numerical computations 
in Sect. IV. The equilibrium of an axisymmetric, self-gravitating perfect fluid is described by the 
gravitational potential density p, pressure p and a 3- velocity 

v" = n^f", (15) 

where 93" is the rotational Killing vector field. 



We will use a Lagrangian perturbation formalism ( [Friedman and Schutz 1978a ) in which 



perturbed quantities are described in terms of a Lagrangian displacement vector that connects 
fluid elements in the equilibrium and perturbed star. The Eulerian change 5Q in a quantity Q is 
related to its Lagrangian change AQ by 

AQ = 5Q + £^Q, (16) 

with £^ the Lie derivative along 

The fluid perturbation is then determined by the displacement 

Av"" = dte (17) 

^ = ^ = -V.r (18) 
IP P 



*Note that for spherical stars, nonhnear coupUngs invahdate the linear approximation after a time t ~ R/5v, 
comparable to the time for a fluid element to move once around the star. For nonzero angular velocity, the linear 
approximation is expected to be valid for all times, if the amplitude is sufficiently small, roughly, if \Sv\ < RQ. 



-7- 



Since the equilibrium spacetime is stationary and axisymmetric, we may decompose our 
perturbations into modes of the form|^ ^i{cTt+rrup) rpj^g corresponding Eulerian changes are 

Sv" = i{a + 771^1)^" (19) 
Sp = -Vaipn (20) 
= 'j^^Sp; (21) 

and the change in the gravitational potential is determined by 

V^J^ = AnGdp. (22) 



We can expand the perturbed fluid velocity, Sv"', in vector spherical harmonics ( |Regge and 
Wheeler 1957, see also p?horne 1980D , 



dv" = l-WiYi'^Vr + ViVYi"" - iUie''''''VbYrVcr\ e'"\ (23) 

and examine the perturbed Euler equation. 

The Lagrangian perturbation of Euler's equation is 

= idt + £v)AVa + Va[A{h-^V^ + <P)], (24) 

and its curl, which expresses the conservation of circulation for an isentropic star, is 

= q'' = i{a + mn)e'''"'VbAvc, (25) 

or 

= q^ = i{a + mn)e''''^VbdVc + Vte'^^^'V b{£ &vV c) ■ (26) 
Using the spherical harmonic expansion (^) of dv"" we can write the components of q"" as 

1 °° r 

= g'' = — ^ [(^^ + mn)l{l + 1) - 2ma]UiY{^ - 2^Vi[sin OdeVr + l{l + 1) cos OYl^] 

^ l=m ^ 

+ 2nWi [sin edeYi"' + 2 cos OYf] | e^'^* , (27) 



^We will always choose m > since the complex conjugate of an m < mode with frequency ct is an m > mode 
with frequency —a. Note that a is the frequency in an inertial frame. 
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= q'' 



J.2 gjj^ 



£ I m(cj + mS7) (drVi - rr - 2ndrVi cos 6 sin BdeYf- 

l=m 



and 



= a^ 



+ 20m2— y,™ - 2ndrWi sin^ (9y,"" - 2mQ.drUi cos 
r 

r 



■ ^ i m{a + m^l)drUiY{^ - l^dr-Ui cos9sm9deYr 

l=m ^ 

+ 2n—[m^ - l{l + 1) sin^ ^]y™ - 2mndrVi cos 9Yf' 
r 



(28) 



+ 



/ Wj\ Vi 



sin 0df,Yr [ e'^K 



(29) 



These components are not independent. The identity VaQ"^ = 0, which fohows from equation ( p5| ) 
serves as a check on the right-hand sides of (p7|) - (|29|). 



(30) 
(31) 



Let us rewrite these equations making use of the standard identities, 



sin^a^y;" 
cos eY{' 



where 



Qi = 



iQi+iYi+i -{1 + mYi-i 

Qi+,Y,^, + QiY,^, 



(l + m){l — m) 



X2l - 1){21 + 1) 
Defining a dimensionless comoving frequency 

(cr + mi}) 



K = 



we find that the = equation becomes 

oo f 

0=E 1 i^i^Kl + 1) - m]UiY- 



l=m 



+{Wi - lVi){l + 2)Q,+iy™, - [Wi + + l)Vi]{l - l)QiYl^^ 



(32) 



(33) 



(34) 



(7 =0 becomes 

cx: 

= E ^ -Qi+iQi+2 



1=771 
+ 



ivl - wl 



^1+2 ~ Ql+1 



(m — hnl)Ui — ml — 



vm 



(i^m + + l)Qf - IQl,) Vl - (l - - gf+i) Wl - \ 



W, 



Km h m 



2V1 



Yr 
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Qi 



(m + ^k{1 + 1)) U[ + m{l + 1 



r 



+QI-1QI 

and q'^ = becomes 



{l + l)V{ + Wl 



(35) 



00 

E 

l=m 



-IQi+iQi- 



U[-{1 + 1 



Ul 



Y, 



+ Q1+1 



[\k.I - m)Vl + ml 



Vi 



1+2 
Wi 



vm 



+ 



(1^^ + (/ + 1)QF - /QF+i) [// + (m^ - /(/ + 1) (1 - q2 _ q2^^ 
+ 1) + m) y/ + m{l + 1)-?^ - + 1)-^ 



r 



Qi 



+ il + l)Qi.iQi 



vm 



r 



vm 



(36) 



where ' = -r- 



dr ' 



From this last form of the equations it is clear that the rotation of the star mixes the axial 
and polar contributions to That is, rotation mixes those terms in ( |23| ) whose limit as — > 
is axial with those terms in ( p3D whose limit as — > is polar. It is also evident that the axial 
contributions to Sv"" with / even mix only with the odd / polar contributions, and that the axial 
contributions with / odd mix only with the even / polar contributions. In addition, we prove in 
appendix A that for non-axisymmetric modes the lowest value of / that appears in the expansion 
of dv"" is always / = m (When m = this lowest value of / is either or 1.) 

Thus, we find two distinct classes of mixed, or hybrid, modes with definite behavior under 
parity. This is to be expected because a rotating star is invariant under parity. Let us call 
a non-axisymmetric0 mode an "axial- led hybrid" (or simply "axial- hybrid" ) if 6v"' receives 
contributions only from 



axial terms with I 
polar terms with / 



m, m + 2, m + 4, ... and 
m + 1, m + 3, m + 5, . . .. 



Such a mode has parity (—1)™+^. 

Similarly, we define a non-axisymmetricQ mode to be a "polar-led hybrid" (or "polar-hybrid" ) 



When m — there exists a set of modes with parity +1 that may be designated as "axial- led hybrids" since Sv"" 
receives contributions only from axial terms with I — 1, 3, 5, . . . and polar terms with Z = 2, 4, 6, . . .. 

^When m — there exist two sets of modes that may be designated as "polar-led hybrids." One set has parity 
— 1 and fJn" receives contributions only from polar terms with Z = 1, 3, 5, . . . and axial terms with Z = 2, 4, 6, . . .. The 
other set (which includes the radial oscillations) has parity and 5v°' receives contributions only from polar terms 
with Z = 0, 2, 4, . . . and axial terms with Z = 1, 3, 5, . . .. 
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if Sv"" receives contributions only from 



polar terms with / = m, m + 2, m + 4, ... and 
axial terms with / = m + 1, m + 3, m + 5, .... 

Such a mode has parity (—1)™'. 

Let us rewrite the equations one last time using the orthogonality relation for spherical 
harmonics, 



mm' 1 



where dD. is the usual solid angle element. 

From equation (|3^ ) we find that / q^Y^'^dVt = gives 
= + 1) - m]Ui + + l)Qi[Wi^i -{I- l)Vi.{\ - IQi+AWi+i + + 2)Vi+i] 



Similarly, / q^Yi*'^dQ. = gives 

= QiQi-i{{l-2)VU-WU} + Qi^[m-\K{l-l)]UU-m{l-l 

+ (l - Q? - Q'+i) Wl - [\Km + + _ /q2^^] vl + l^rn^ - rn^}l 
+ Qi+i I [m + ^k{1 + 2)]C//+i + m{l + 2)^ 

-Qi+2Qi+i{ii + m'+2 + w[+2} 

and / q'^Yi*"^dn = gives 



-{l-2)QiQ 



i-i 



U, 



1-2 



(l-l) 



u, 



1-2 



+ 



+ Ql \ [h<l - 1) - m]Vi'_i + m{l - 1) 



+ {l + 3)Qi+2Qi+i 

- l{l + 1) (l - Q? - Qf+i 



Ul+2 



r 



Qi+i + 2) + m]y4i + mil + 2)^ - + 2)^| . 



(37) 



(38) 



(39) 



(40) 



4. Method of Solution 

In our numerical solution, we restrict consideration to slowly rotating stars, finding axial- and 
polar-led hybrids to lowest order in the angular velocity 17. That is, we assume that perturbed 
quantities introduced above obey the following ordering in 0: 



Wi ~ 0(1), Vi - 0(1), Ui - 0(1), 
5/9~0(17), 5p~0(S7), 5«>~0(0), cj~0(Jl). 



(41) 
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The ^ limit of such a perturbation is a sum of the zero- frequency axial and polar perturbations 
considered in Sect. II. Note that, although the relative orders of 6p and (5t>° are physically 
meaningful, there is an arbitrariness in their absolute order. If (Sp,6v"') is a solution to the 
linearized equations, so is (Q5p,Q5v'^). We have chosen the order ( |4l| ) to reflect the existence 
of well-defined, nontrivial velocity perturbations of the spherical model. Other authors (e.g., 
Lindblom and Ipser (1998)) adopt a convention in which 5v^ = 0(i7) and 5p = 0{Q'^). 

To lowest order, the equations governing these perturbations are the perturbed Euler 
equations (38) - (|40|) and the perturbed mass conservation equation, (^, which becomes 



rWl + (^1 + r^"^ Wi - l{l + l)Vi = 0. (42) 

In addition, the perturbations must satisfy the boundary conditions of regularity at the center 
of the star, r = 0, regularity at the surface of the star, r = R, and the vanishing of the Lagrangian 
change in the pressure at the surface of the star, 

= Ap = 6p + £i:p = C''p' + O{Q). (43) 



Equations (|38| ) through (42) are a system of ordinary differential equations for Wi'{r), Vii{r) 
and Ui'{r) (for all /'). Together with the boundary conditions, these equations form a non-linear 
eigenvalue problem for the parameter k, where /tO is the mode frequency in the rotating frame. 

To solve for the eigenvalues we proceed as follows. We first ensure that the boundary 
conditions are automatically satisfied by expanding Wi'{r), Vi'{r) and Ui'{r) (for all /') in regular 
power series about the surface and center of the star. Substituting these series into the differential 
equations results in a set of algebraic equations for the expansion coefficients. These algebraic 
equations may be solved for arbitrary values of k using standard matrix inversion methods. For 
arbitrary values of k, however, the series solutions about the center of the star will not necessarily 
agree with those about the surface of the star. The requirement that the series agree at some 
matching point, < vq < R, then becomes the condition that restricts the possible values of the 
eigenvalue, kq. 

The equilibrium solution (/?, <I>) appears in the perturbation equations only through the 
quantity {p' / p) in equation (^). We begin by writing the series expansion for this quantity about 
r = as 

i odd 

and about r = R as 

where the tTj and TTk are determined from the equilibrium solution. 
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Because ( pq ) relates Ui{r) algebraically to Wi±i{r) and Vi±i{r), we may eliminate Uii{r) (all 
I') from (|3^) and ([40|). We then need only work with one of equations ( p9| ) or ( ^0[ ) since the 
equations (38) through (EH) are related by VaQ"' = 0. 



We next replace p'/p, Wi', and Vi' in equations ( pOD or ( [40| ) by their series expansions. We 
eliminate the Ui'{r) from either ( p9| ) or ( |40D and, again, substitute for the W/'(r) and Vi/{r). 
Finally, we write down the matching condition at the point ro equating the series expansions about 
r = to the series expansions about r = R. (Explicitly one equates ( p6[) and ( p7[) of Appendix 



B for axial-led modes or ( B13 ) and ( B14 ) for polar-led modes). The result is a linear algebraic 



system which we may represent schematically as 

Ax = 0. (46) 

In this equation, ^ is a matrix which depends non-linearly on the parameter k, and x is a vector 
whose components are the unknown coefficients in the series expansions for the W/'(r) and Vii{r). 
In Appendix B, we explicitly present the equations making up this algebraic system as well as the 
forms of the regular series expansions for VF/'(r) and VJ/(r). 

To satisfy equation (^) we must find those values of k for which the matrix A is singular, i.e., 
we must find the zeroes of the determinant of A. We truncate the spherical harmonic expansion 
of dv"" at some maximum index l,^^^ and we truncate the radial series expansions about r = and 
r = ii at some maximum powers z,„ax and k^^^, respectively. 

The resulting finite matrix is band diagonal. To find the zeroes of its determinant we use 
standard root finding techniques combined with routines from the LAPACK linear algebra libraries 
( Anderson et al. 19941) . We find that the eigenvalues, kq, computed in this manner converge 



quickly as we increase Z^axi Vax and k,^^^. 

The eigenfunctions associated with these eigenvalues are determined by the perturbation 
equations only up to normalization. Given a particular eigenvalue, we find its eigenfunction by 
replacing one of the equations in the system (^) with the normalization condition that 

ymii" = R) = 1 for polar-hybrids, or that ^^^^ 
Vm+iir = R) = 1 for axial- hybrids. 

Since we have eliminated one of the rows of the singular matrix A in favor of this condition, the 
result is an algebraic system of the form 

Ax = b, (48) 

where A is now a non-singular matrix and 6 is a known column vector. We solve this system for 
the vector x using routines from LAPACK and reconstruct the various series expansions from this 
solution vector of coefficients. 
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5. The Eigenvalues and Eigenfunctions 

We have computed the eigenvalues and eigenfunctions for uniform density stars and for n = 1 
polytropes, models obeying the polytropic equation of state p = Kp^, where K is a constant. 
Our numerical solutions for the uniform density star agree with the recent results of Lindblom 
and Ipser (1998) who find analytic solutions for the hybrid modes in rigidly rotating uniform 
density stars with arbitrary angular velocity - the Maclaurin spheroids. Their calculation uses the 
two-potential formalism ( Ipser and Managan 1985 ; and Ipser and Lindblom 1990| ) in which the 



equations for the perturbation modes are reformulated as coupled differential equations for a fluid 
potential, 6U, and the gravitational potential, (5$. All of the perturbed fluid variables may be 
expressed in terms of these two potentials. The analysis follows that of Bryan (1889) who found 
that the equations are separable in a non-standard spheroidal coordinate system. 

The Bryan/Lindblom-Ipser eigenfunctions 6Uq and turn out to be products of associated 
Legendre polynomials of their coordinates. This simple form of their solutions leads us to expect 
that our series solutions might also have a simple form - even though their unusual spheroidal 
coordinates are rather complicated functions of r and 0. In fact, we do find that the modes of 
the uniform density star have a particularly simple structure. For any particular mode, both the 
angular and radial series expansions terminate at some finite indices Iq and io (or ko)- That is, the 
spherical harmonic expansion (^) of dv"" contains only terms with m < I < Iq for this mode, and 
the coefficients of this expansion - the Wi{r), Vi{r) and Ui{r) - are polynomials of order Iq. For all 
Iq > m there exist a number of modes terminating at Iq. 

In Tables || to ^ we present the functions Wi{r), Vi{r) and Ui{r) for all of the axial- and 
polar- led hybrids with m = 1 and m = 2 for a range of values of the terminating index Iq. (See 
also Figure [l|.) For given values of m > and Iq there are Iq — m + 1 modes. (When m = there 
are rnodes. See equation (^) below.) We also find that the last term in the expansion (pSj), the 
term with / = Iq, is always axial for both types of hybrid modes. This fact, together with the fact 
that the parity of the modes is, 

(—1)™ for polar- led hybrids 
(—1)™+-^ for axial-led hybrids, 

(for m > 0) implies that Iq — m+\ must be even for polar-led modes and odd for axial-led modes. 



The fact that the various series terminate at Iq, iq and /cq implies that Equations ([46|) and 
will be exact as long as we truncate the series at /^ax > ^o, ^max > io and A;,nax > ^o- 

To find the eigenvalues of these modes we search the k axis for all of the zeroes of the 
determinant of the matrix A in equation (|4^). We begin by fixing m and performing the search 
with l^^^ = m. We then increase /,„ax by 1 and repeat the search (and so on). At any given value 
of ^max) the search finds all of the eigenvalues associated with the eigenfunctions terminating at 

^0 — ^max* 

In Table ^, we present the eigenvalues kq found by this method for the axial- and polar-led 
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hybrid modes of uniform density stars for a range of values of Iq and m. Observe that many 
of the eigenvalues, (marked with a *) satisfy the condition a{a + mfi) < 0. [Recall that the 
mode frequency in an inertial frame is o" = {kq — m)0,]. The modes whose frequencies satisfy this 
condition are subject to a gravitational radiation driven instability in the absence of viscosity. The 
modes having Iq = m > (or Iq = 1 for m = 0) are the purely axial r-modes. Their frequencies 
were found by Papalouizou and Pringle (1978) and are given by kq = 2/{m + 1) (or hq = for 
m = 0). We find that there are no purely polar modes satisfying our assumptions (|4l| ) in these 
stellar models. 

We have compared these eigenvalues with those of Lindblom and Ipser (1998). To lowest 
non-trivial order in Q their equation for the eigenvalue, kq, can be expressed in terms of associated 
Legendre polynomials^ (see Lindblom and Ipser's equation 6.4), as 

(4 - ^2)^p-Yi(f ) - 2m/^™ i(f ) = 0. (50) 

For given values of m > and Iq this equation has Iq — m -\- 1 roots (corresponding to the number 
of distinct modes), which can easily be found numerically. (For m = there are Iq roots.) For 
the range of values of m and Iq checked our eigenvalues agree with these to machine precision. 
(Compare our Table |5| with Table 1 in Lindblom and Ipser 199^ .) 



We have also compared our eigenfunctions with those of Lindblom and Ipser. For a uniformly 
rotating, isentropic star, the fluid velocity perturbation, 6v"', is related (Ipser and Lindblom 1990) 
to 5U by 

VJU = - [iKQgab + 2VbVa] 5v\ (51) 
Since the </? component of this equation is simply 

"2 



imSU = —Q.r^ sin^ 



■5v'' + 2 cot eSv'^ + iK5v'^ 
r 



(52) 



it is straightforward numerically to construct this quantity from the components of our dv"" and 
compare it with the analytic solutions for 5U given by Lindblom and Ipser (see their equation 
7.2). We have compared these solutions on a 20 x 40 grid in the (r — 9) plane and found that they 
agree (up to normalization) to better than 1 part in 10^ for all cases checked. 

Because of the use of the two-potential formalism and the unusual coordinate system used 
in their analysis, the axial or polar hybrid character of the Bryan/Lindblom-Ipser solutions is 
not obvious. Nor is it evident that these solutions have, as their $7 — > limit, the zero-frequency 
convective modes described in Sect. II. The comparison of their analytic results with our numerical 
work has served the dual purpose of clarifying these properties of the solutions and of testing 
the accuracy of our code. The computational differences are minor between the uniform density 



The index / used by Lindblom and Ipser is related to our Iq by Z = /q + 1- Our convention agrees with the usual 
labelling of the Iq = m pure axial modes. 
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calculation and one in which the star obeys a more realistic equation of state. Thus, this testing 
gives us confidence in the validity of our code for the polytrope calculation. As a further check, 
we have written two independent codes and compared the eigenvalues computed from each. One 
of these codes is based on the set of equations described in Appendix B. The other is based on 
the set of second order equations that results from using the mass conservation equation, (^), to 
substitute for all the Vi{r) in favor of the Wi{r). 

For the n = 1 polytrope we will consider and, more generally, for any isentropic equation 
of state, the purely axial r-modes are independent of the equation of state. In both isentropic 
and non-isentropic stars, pure r-modes exist whose velocity field is, to lowest order in 0,, 
an axial vector field belonging to a single angular harmonic (and restricted to harmonics 
with / = m in the isentropic case). The frequency of such a mode is given (to order by 
= {a + mQ) = 2mQ/l{l + 1) ( Papalouizou and Pringle 1978D and is independent of the equation 



of state. In isentropic stars, only those modes having / = m (or / = 1 for m = 0) exist, and for 
these modes the eigenfunctions are also independent of the (isentropic) equation of state|^. This 
independence of the equation of state occurs for the r-modes because (to lowest order in Q) fluid 
elements move in surfaces of constant r (and thus in surfaces of constant density and pressure). 
For the hybrid modes, however, fluid elements are not confined to surfaces of constant r and one 
would expect the eigenfrequencies and eigenfunctions to depend on the equation of state. 

Indeed, we find such a dependence. The hybrid modes of the n = 1 polytrope are not 
identical to those of the uniform density star. On the other hand, the modes do not appear to 
be very sensitive to the equation of state. We have found that the character of the polytropic 
modes is similar to the modes of the uniform density star, except that the radial and angular 
series expansions do not terminate. For each eigenfunction in the uniform density star there is a 
corresponding eigenfunction in the polytrope with a slightly different eigenfrequency (See Table |6|.) 
For a given mode of the uniform density star, the series expansion (|2^) terminates at I = Iq. For 
the corresponding polytrope mode, the expansion ( [2^ ) does not terminate, but it does converge 
quickly. The largest terms in ( p3| ) with / > Iq are more than an order of magnitude smaller than 
those with / < Iq and they decrease rapidly as I increases. Thus, the terms that dominate the 
polytrope eigenfunctions are those that correspond to the non-zero terms in the corresponding 
uniform density eigenfunctions. 

In Figures |l] and |2| we display the coefficients Wi{r), V/(r) and Ui{r) of the expansion (p3|) for 
the same m = 2 axial-led hybrid mode in each stellar model. For the uniform density star (Figure 
I) the only non-zero coefficients for this mode are those with ^ < = 4. These coefficients are 
presented explicitly in Table ^ and are low order polynomials in r. For the corresponding mode 
in the polytrope, we present in Figure I the first seven coefficients of the expansion (Eq). Observe 



et al 



®The only non-zero term in ( p3[ ) for these modes is the axial lo = m term with coefficient Um{r) — r"^'^^ (Provost 
f98f ). The purely axial mode with m = has lo = 1 and radial dependence r^. This mode has zero frequency 



and corresponds to a small uniform change in the angular velocity of the star. 
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that those coefficients with I < 4 are similar to the corresponding functions in the uniform density 
mode and dominate the polytrope eigenfunction. The coefficients with 4 < I < 6 are an order of 
magnitude smaller than the dominant coefficients and those with I > 6 are smaller still. (Since 
they would be indistinguishable from the {r/R) axis, we do not display the coefficients having 
/ > 6 for this mode.) 

Just as the angular series expansion fails to terminate for the polytrope modes, so too do 
the radial series expansions for the functions Wi{r), ^lif) and Ui{r). We have seen that in the 
uniform density star these functions are polynomials in r (Tables Q through In the polytropic 
star, the radial series do not terminate and we are required to work with both sets of radial series 
expansions - those about the center of the star and those about its surface - in order to represent 
the functions accurately everywhere inside the star. 

In Figures |^ through 11 we compare corresponding functions from each type of star. For 
example. Figures ^, ^, and ^ show the functions Wi{r), Vi{r) and f//(r) (respectively) for / < 6 for 
a particular m = 1 polar-led hybrid mode. In the uniform density star this mode has eigenvalue 
kq = 1.509941, and in the polytrope it has eigenvalue kq = 1-412999. The only non-zero functions 
in the uniform density mode are those with I <Iq = 2 and they are simple polynomials in r (see 
Table |2|). Observe that these functions are similar, but not identical to, their counterparts in the 
polytrope mode, which have been constructed from their radial series expansions about r = and 
r = R (with matching point tq = 0.5R). Again, note the convergence with increasing / of the 
polytrope eigenfunction. The mode is dominated by the terms with / < 2 and those with / > 2 
decrease rapidly with /. (The / = 5 and / = 6 coefficients are virtually indistinguishable from the 
{r/R) axis.) 

Because the polytrope eigenfunctions are dominated by their / < Iq terms, the eigenvalue 
search with /^^x = will find the associated eigenvalues approximately. We compute these 
approximate eigenvalues of the polytrope modes using the same search technique as for the 
uniform density star. We then increase /^ax and search near one of the approximate eigenvalues for 
a corrected value, iterating this procedure until the eigenvalue converges to the desired accuracy. 
We present the eigenvalues found by this method in Table |6|. 

As a further comparison between the mode eigenvalues in the polytropic star and those in 
the uniform density star we have modelled a sequence of "intermediate" stars. By multiplying the 
expansions (U) and (^) for {p' / p) by a scaling factor, e G [0,1], we can simulate a continuous 
sequence of stellar models connecting the uniform density star (e = 0) to the polytrope (e = 1). 
We find that an eigenvalue in the uniform density star varies smoothly as function of e to the 
corresponding eigenvalue in the polytrope. 
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6. The Effects of Dissipation 



The effects of gravitational radiation and viscosity on the pure Iq = m r-modes have already 
been studied by a number of authors. (Lindblom et al 1998, Pwen et al. 1998 , Andersson et al 
1998, Kokkotas and Stergioulas 1998| , [Lindblom et al. 1999| ) All of these modes are unstable to 
gravational radiation reaction, and for some of them this instability strongly dominates viscous 
damping. We now consider the effects of dissipation on the axial- and polar-hybrid modes. 

To estimate the timescales associated with viscous damping and gravitational radiation 



reaction we follow the methods used for the Iq = m modes (Lindblom et al 1998, see also [pser 
and Lindblom 1991| ) . When the energy radiated per cycle is small compared to the energy of the 
mode, the imaginary part of the mode frequency is accurately approximated by the expression 



1 
r 



1 dE 
'2E~dt' 



(53) 



where E is the energy of the mode as measured in the rotating frame, 

' 5p 



E 



cPx. 



(54) 



The rate of change of this energy due to dissipation by viscosity and gravitational radiation is. 



dE 



l>2 



(55) 



The first term in ( |55D represents dissipation due to shear viscosity, where the shear, 6aab, of the 
perturbation is 

S(^ab = I (Va^Vfe + VbSVa - lOab^cSv^) , (56) 

and the coefficient of shear viscosity for hot neutron-star matter is ( Cutler and Lindblom 1987] ; 



Sawyer 19891) 



2 X 10 
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P 
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— Tf, — g-cm -s . 



(57) 



The second term in ( pSD represents dissipation due to bulk viscosity, where the expansion, 66, 
of the perturbation is 



69 = Vr6v' 



(58) 



and the bulk viscosity coefficient for hot neutron star matter is (Cutler and Lindblom 1987; Sawyer 

6 



1989) 



C = 6 X 10 



25 
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a + mnj VlO^^g-cm-V \W^K 
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g-cm ^-s ^. 



(59) 
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The third term in ( |55| ) represents dissipation due to gravitational radiation, with couphng 
constant 

4^G (/ + !)(/ + 2) 



Ni 



r,2l + l 



l{l - 1)[{21 + IW 



(60) 



The mass, SDim, and current, 6 Jim, multipole moments of the perturbation are given by ( p?horng 

(61) 



1980, Lindblom et al 1998) 



and 



c\l + l 



r' {p6Va + 6pva) Yi''f*d^x 



where is the magnetic type vector spherical harmonic ( Thorne 1980| ) given by. 



Y, 



a,B 



Im 



(62) 



(63) 



To lowest order in 0, the energy ( ^4|) of the hybrid modes is positive definite. Their stability 
is therefore determined by the sign of the right hand side of equation (|55|). We have seen that 
many of the hybrid modes have frequencies satisfying a{a + mJ7) < 0. This makes the third 
term in (^) positive, implying that gravitational radiation reaction tends to drive these modes 
unstable. ( |Chandrasekhar 1970 ; Friedman and Schutz 1978b ; Friedman 197^ ) To determine the 
actual stability of these modes, we must evaluate the various dissipative terms in (|55[). 



We first substitute for 5v'^ the spherical harmonic expansion (23) and use the orthogonality 
relations for vector spherical harmonics ( [Thorne 198C| ) to perform the angular integrals. The 
energy of the modes in the rotating frame then becomes 



E 



oo 

Eh 

l=m 



R 



wf+i{i + i)v;' + i{i + i)ui 



dr. 



(64) 



To calculate the dissipation due to gravitational radiation reaction we must evaluate the 
multipole moments (|6^) and (|62|) . To lowest order in Q the mass multipole moments vanish and 
the current multipole moments are given by 



A 7 2^ 

OJlm — — 



R 



(65) 



To calculate the dissipation due to bulk viscosity we must evaluate the expansion, 56 = S/cSv'^, 
of the perturbation. For uniform density stars this quantity vanishes identically by the mass 
conservation equation For the Iq = m, pure axial modes the expansion, again, vanishes 
identically, regardless of the equation of state. To compute the bulk viscosity of these modes it is 
necessary to work to higher order in (Andersson et al. 1998, Lindblom et al. 1999). On the other 
hand, for the new hybrid modes in which we are interested, the expansion of the fluid perturbation 
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is non-zero in the slowly rotating polytropic stars. After substituting for Sv"^ its series expansion 
and performing the angular integrals, the bulk viscosity contribution to ([55|) becomes 



dE\ 



B 



(66) 



In a similar manner, the contribution to (^) from shear viscosity becomes 



dE 



dtl- ^ L ^2 



l=m 



R 



27] 







+ 1){2P + 21- 3)V;2 + l{l + l)Wi 



+ + +l-2)Uf 




Given a numerical solution for one of the hybrid mode eigenfunctions, these radial integrals 
can be performed numerically. The resulting contributions to ( [55| ) also depend on the angular 
velocity and temperature of the star. Let us express the imaginary part of the hybrid mode 
frequency (^) as. 



1 

r 




(68) 



where p is average density of the star. (Compare this expression to the corresponding expression 
in Lindblom et al. (1998) - their equation (22) - for the Iq = m pure axial modes.) 

The bulk viscosity term in this equation is stronger by a factor than that for the Iq = m 
pure axial modes. This is because the expansion 66 of the hybrid mode is nonzero to lowest order 
in Q for the polytropic star, whereas it is order 0^ for the pure axial modes. This implies that the 
damping due to bulk viscosity will be much stronger for the hybrid modes than for the pure axial 
modes in slowly rotating stars. 



Note that the contribution to ( p8[) from gravitational radiation reaction consists of a sum over 
all the values of / with a non-vanishing current multipole. This sum is, of course, dominated by 
the lowest contributing multipole. 

In Tables to we present the timescales for these various dissipative effects in the uniform 
density and polytropic stellar models that we have been considering with R = 12.57km and 
M = IAMq. For the reasons discussed above, we do not present bulk viscosity timescales for the 
uniform density star. 

Given the form of their eigenfunctions, it seems reasonable to expect that some of the unstable 
hybrid modes might grow on a timescale which is comparable to that of the pure Iq = m r-modes. 
For example, the m = 2 axial-led hybrids all have U2{r) ^ (see, for example. Figures |^ and ^). 
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By equation (65), this leads one to expect a non-zero current quadrupole moment 5J22, and this 
is the multipole moment that dominates the gravitational radiation in the r-modes. Upon closer 
inspection, however, one finds that this is not the case. In fact, we find that all of the multipoles 
5 Jim vanish (or nearly vanish) for I < Iq, where Iq is the largest value of I contributing a dominant 
term to the expansion (23) of 5v^. 



In the uniform density star, these multipoles vanish identically. Consider, for example, the 
m = 2, Iq = A axial-hybrid with eigenvalue k = 0.466901. (See Table ^ and Figure ^ For this 
mode, U2 oc {7x^ — 9x^), where x = (r/R). By equation (|65|), we then find that 

SJ22 oc / x^{7x^ - 9x^)dx = 0, (69) 
Jo 

and that 5J42 is the only non-zero radiation multipole. In general, the only non-zero multipole for 
an axial- or polar- hybrid mode in the uniform density star is dJi^m- 

That this should be the case is not obvious from the form of our eigenfunctions. However, 
Lindblom and Ipser's (1998) analytic solutions provide an explanation. Their equations (7.1) and 
(7.3) reveal that the perturbed gravitational potential, 5$, is a pure spherical harmonic to lowest 
order in 0,. In particular, 

5c^ocy,-i. (70) 

This implies that the only non-zero current multipole is SJi^m- 

We find a similar result for the polytropic star. Because of the similarity between the modes in 
the polytrope and the modes in the uniform density star, we find that although the lower I current 
multipoles do not vanish identically, they very nearly vanish and the radiation is dominated by 
higher / multipoles. 

The fastest growth times we find in the hybrid modes are of order 10^ seconds (at 10^-fC 
and Vl = yJnGp). Thus, the spin-down of a newly formed neutron star will be dominated by the 
Iq = m = 2 mode with contributions from the Iq = m pure axial modes with 2 < m ^ 10 and from 
the fastest growing hybrid modes. 



7. Discussion 

There is substantial uncertainty in the cooling rate of neutron stars, with rapid cooling 
expected if stars have a quark interior or core, or a kaon or pion condensate. Madsen (1998) 
suggests that an observation of a young neutron star with a rotation period below 5 — 10ms would 
be evidence for a quark interior; but even without rapid cooling, the uncertainty in the superfluid 
transition temperature would allow a superfluid to form at about IQ^^K, killing the instability. 
The nonaxisymmetric instability has been expected not to play a role in old neutron stars spun 
up by accretion, because of the high shear viscosity associated with an expected temperature 
< lO^K; but even this is not certain ( [Andersson, Kokkotas and Stergioulas 1998| ). 
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An extension of our numerical method to find modes of rapidly rotating Newtonian models 
and slowly rotating relativistic models appears feasible. Work is in progress to understand the 
way in which the modes join the r- and g- modes of stars that are not isentropic ( [Andersson et al" 
1999D. 
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A. The character of the modes of rotating isentropic stars 

For an equilibrium model that is axisymmetric and invariant under parity, one can resolve 
any degeneracy in the perturbation spectrum to make each discrete mode an eigenstate of parity 
with angular dependence e*™*^. The following theorem holds. 

Theorem 1 Let {dp^Sv"") with dv"" ^ be a discrete normal mode of a uniformly rotating stellar 
model obeying a one-parameter equation of state. Then the decomposition of the mode into 
spherical harmonics Yj^ (i.e., into {l,m) representations of the rotation group about its center of 
mass) has I = m as the lowest contributing value of I, when m ^ 0; and has or 1 as the lowest 
contributing value of I, when m = 0. 

In Sect. Ill, we designate non-axisymmetric modes with parity (—I)™' "polar-led hybrids", 
and non-axisymmetric modes with parity (—1)™''"^ "axial-led hybrids," and briefly discuss the 
m = case. 

Note that the theorem holds for p-modes as well as for the rotational modes that are our 
main concern. A p-mode is determined by its density perturbation and is therefore dominantly 
polar in character regardless of its parity. For a rotational mode, however, the lowest / term in its 
velocity perturbation is at least comparable in magnitude to the other contributing terms. 

We prove the theorem separately for each parity class. 

A.l. Axial-Led Hybrids with m > 

Let I be the smallest value of for which Ui' 7^ in the spherical harmonic expansion (^) of 
the perturbed velocity field The axial parity of (—1)'"''^, and the vanishing of YJ™ for 
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I < m implies / > m. That the mode is axial-led means Wii = and Vii = for /' < I. We show by 
contradiction that I = m. 

Suppose I > m + I. From equation (|3^), / q^Y^^dQ = 0, we have 

+ 1) - m]Ui = lQi+i[Wi+i + (/ + 2)Vz+i], (Al) 
and from equation ( ^9|) with / replaced by / — 1, / q^Yj*"^dn = 0, we have 

Qi+i [(/ + 2)y/+i + W,Vi] = {[m + + !)][// + m(/ + 1)^1 . (A2) 
These two equations, together imply that 

U[ + -Ui = 0, 
r 

or 

C/i = Kr-\ 

which is singular at r = 0. 



A. 2. Axial-Led Hybrids with m = 

Let m = and let I be the smallest value of /' for which C/p 7^ in the spherical harmonic 
expansion ( |2^ ) of the perturbed velocity field dv"". Since VaYQ = 0, the mode vanishes unless 
/ > 1. That the mode is axial- led means Wi' = and Vi' = for V <l. We show by contradiction 
that I = 1. 

Suppose / > 2. Then / q^Yi^d^ = becomes, 

Ul + -Ui = Q, (A3) 
r 

or 

Ui = Kr-\ 

which is singular at r = 0. 



A. 3. Polar-Led Hybrids with m > 

Let / be the smallest value of /' for which Wii 7^ or V// 7^ in the spherical harmonic 
expansion (p3|) of the perturbed velocity field Sv"" . The polar parity of Sv*^ , {~^)\ and the 
vanishing of Yj^ for / < m implies / > m. That the mode is polar-led means Ui' = for I' < I. We 
show by contradiction that I = m. 
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Suppose I > m + 1. Then / q^Y^^dD. = becomes 

Wi + {1 + l)Vi = 0, 

and / qfYi*!^dn = becomes, 

= - + 1) + m]V{ + m{l + 1)^ - ^k{1 + l)^} 



These two equations, together imply that 



1) +m] 



^/ + (/ + !)- 



+ (/ + 2)Q/+i 



0, 



or 



a^il + 1) + rn]Vi + + 2)Qi+iUi+i = Kr-^^+^\ 



(A4) 



(A5) 



which is singular at r = 0. 



A. 4. Polar-Led Hybrids with m = 

Let m = and let I be the smallest value of /' for which Wii ^ and VJ/ 7^ in the spherical 
harmonic expansion (23) of the perturbed velocity field 5f When / = the mode is automatically 
polar- led; thus we need only consider the case / > 1. That the mode is polar-led means Ui' = for 
I' < I. We show by contradiction that 1 = 1. 



Suppose I > 2. Then / q'^Yi*_^di} = becomes 



Wi + {1 + l)Vi = 0, 



and / g'^yjlV^^ = becomes. 



W, 



+ {1 + 2)Qi+i 



Ui+i 



These two equations, together imply that 



Vi' + il + l 



Vi 



+ (/ + 2)Q;+i 



uUi + {1 + 1 



or 



+ l)Vi + {l + 2)Qi+iUi+i = Kr 



(A6) 



(A7) 



which is singular at r = 0. 
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B. The algebraic equations governing the hybrid modes to lowest order in Q. 

In this appendix, we make use of the following definitions: 



For reference, we repeat the definitions ( |32|) and 



ai = ^Km + {I + l)Qf - IQf^^ (Bl) 
bi ^ rri^ - /(/ + 1) (l - q2 _ g2^^^ (B2) 

Q = ^kI{1 + 1) - m (B3) 



(B4) 



(B5) 



Qi 



n 

[l + m){l — m) 



{21 - 1)(2/ + 1) 



B.l. Axial Hybrids 

For / = m, m + 2, m + 4, . . . the regular series expansions]^ about the center of the star, r = 0, are 

m+j oo 



/ r \ "i+J ^ fry 
Wm+j+i{r) = i^-j 22 "^j+^^i [-^J (B6a: 

i odd 

/ ^ \ m+j oo / r \ ' 



i odd 
i odd 

where j = 0, 2, 4, . . .. 

The regular series expansions about r = R, which satisfy the boundary condition Ap = are 

oo / ^ \ fc 



W^+j+i{r) = J2wj+i,k[l-j,) (B7a) 



k 



y„+j+i(r) = J2 ^i+i,fc ( 1 - (B7b) 



fc=0 



oo / ^ N fc 

[/„+,(r) = ^n,-,(l--j (B7c) 



fc=0 



^"We present the form of the series expansions for Ui{r) for reference; however, we do not need these series since 
we eliminate the Ui{r) using equation (^sj). 
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where j = 0, 2, 4, . . .. 

These series expansions must agree in the interior of the star. We impose the matching condition 



that the series (B6) truncated at i^ax be equal at the point r = ro to the corresponding series (B7) 
truncated at fcmax- That is, 
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i odd 
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i odd 
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^max / \ k 

fc=l V -TL. 



t "-max 



fc=0 



(B8a) 
(B8b) 



When we substitute ( [B6| ) and (|44| ) into (|42|), the coefficient of (r/i?)™^-'^* in the resulting equation 



IS 



i-2 



= (m + j + i + l)wj+i^i + ^ VTs-Wj+i^i^s-i - (m + j + l)(m + j + 2)?jj+i,i (B9) 



s = l 
s odd 



Similarly, when we substitute (|B7| ) and ( ^5| ) into (42), the coefficient of [1 — {r /R)]^ in the resulting 
equation is 

k 

= {k+l) [wj+i^k - Wj+i^k+i\+Y ~ *s-2) Wj+i^k-s+1 - {m + i + 2){m + j + l)vj+i^k (BIO) 

s=Q 

where we have defined 7t_2 = = Wj+i,o- 

When we use (^) to eliminate the Ui{r) from ( |40| ) and then substitute for the Ty/±i(r) and V/±i(r) 
using (p6D, the coefficient of (r/i?)™+-'+* in the resulting equation is 



= {i + l){m + j -2){m+ j -l)Qm+jQm+j-lQm+j-2Cm+jCm+j+2 

X i(;j_3,i+4 - (m + j - 3)t;j_3,i+4 

- (5m+jCm+i+2| (« + l)(m + j - 2fQ'^^j^j_^Cm+j + ^K{m + j - l)Cm+j-2Cm+j 

+ (m + j + 1) [{m + j + i)ara+j + hm+j] Cm+j-2\wj^i^i+2 
+ Qm+jCm+j+2^ 



Cm+j — 2 C-m+j 



^K{m + j — l){m + j + i) — (i + l)m 
+ {m + j + l){m + j - 1) [{m + j + i)am+j + bm+j] Cm+j-2 
- ii + l){m+ j){m+ j - 2fQl^_^j_^Cm+j f^^i-i,j+2 



+ Qm+j+lCrn+j ^2^^1^(771 + j + 2) 

Crn+jCm+j+2 + + j) [{n^ + J ' + *)Or?i+j + 

- (2m + 2j + i + 2)(m + j + 3)2Q^+j+2Cm+i 



Cm+j+2 



+ (5m+j+iCm+j_2< {m+j + 2){m+j) [{m + j + z)a 



|K(m + j + 2)(m + J + i) + m(2m + 2j + i + 2) 



+ (2m + 2j + i + 2)(m + i + 3f{m + j + l)Q^+j-+2Cm+j 
+ (2m + 2i + i + 2)(m + i + 3)(m + j + 2)Qm+i+3Qm+j+2Q m+j + 1 Cm+j — 2 C-m+j 

X If j+3,i_2 + (m + j + 4)t;j+3_j_2 



(Bll) 



When we use (^) to eliminate the [/^(r) from ( |40| ) and then substitute for the Ty;±i(r) and V/±i(r) 
using (pTj), the coefficient of [1 — (r/i?)]'^ in the resulting equation is 



= -(m + j - A; - l)(m + j - l)(m + j - 2)Q„i+jQrn+j-lQm+j~-2Cm+jCm+j+2 

Wj-3,k - {m + j - 3)vj-3^k 



X 



- {k + l)(m + j - l)(m + j - 2)Qm+jQm+j-lQm+j-2Cm+jCm+j+2 

X Wj-3^k+i - {m + j - 3)5j_3,fc+i 

+ Qm+jCm+j+2^^{m + j - k - l){m + j - 2fQ'^_^_j_^Cm+j - \n{m + j - l)Cm+j-2Cm+j 
- {m + j + 1) {hm+j + kam+j) Cm+j-2 \wj-l,k 



+ {k + l)Qm+jCm+j+2 \ {m + j-2fQl+^ -iCm+j + {m + j + l)am+jCm+j~2 >'Wj-l,fc+l 



+ Qm+jCm+j+2\ {m + j - k-l){m + j - 2f{m + i)(5m+j_iCm+j 



+ 



\i^k{m + j — 1) + m{m + j — k — 1) 



Cm+j — 2 Cm+j 



+ {m + j + l){m + j - 1) {brn+j + kam+j) Cm+j-2 \vj-l,k 



+ {k + l)Qm+jCm+j+2 I {m + j){m + j-2YQi^^^ 

— 1 Cm+j 



+ 



m — ^K{m + j — 1) 



C-m+j — 2 Cm+j 
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- (m + j + l)(m + j - l)am+jCrn+j-2 
Qm+j+lCm+j~2^{'m + j) (bjn+j + kam+j) Cm+i+2 + \n{m + j + 2)Cm+jCm+j+2 

- {m + j + k + 2){m + j + 3fQl^^j^2'^rn+j \wj+i,k 



+ l)Qm+j+lCm+j-2'N -{m + j)am+jCrn+j+2 + {m + j + 3)^(5^+j+2Cm+j \'Wj+l,k+l 



+ <3m+j+lCm+j-2< (m + j + 2){m + j) (6m+j + kam+j) Cm+j+2 



m{m + j + k + 2) + ^Kk{m + j + 2) 



+ {m + j + k + 2){m + j + 3)2 (m + j + l)g^+,.+2Q 



Cm+j Cm+j+2 

+j+2^m+j }Vjj^i^k 



+ {k + l)Qm+j+lCm+j-2{ -{m+j + 2){m + j)am+jCm+j+2 



+ 



^K{m + j + 2) + m 



C-m+j Cm+j+2 



im+j+2Cm+j jVj+l^k+l 
+ (m + j + k + 2){m + j + 3){m + j + 2)Qm+j+zQm+j+2Qm+j+lCm+j~2Cm+j 

X *j+3,fc + (m + J + 4)t;j+3,fc 

- (/C + l)(m + j + 3)(m + j + 2)(5m+j+3(5m+j+2(5m+j+lCm+j_2Cm+i 

w;-,-+3,fc+i + (m + j + 4)5j+3^fc+i 



(B12) 



The equations ( p8D through ( |B12| ) make up the algebraic system (46) for eigenvalues of the 
axial-led hybrid modes. One truncates the angular and radial series expansions at indices jmaxj 
iniax and /c^ax and constructs the matrix A by keeping the appropriate number of equations for 
the number of unknown coefficients Wj+i^j, "Wj+i.fc and Vj+i^k- In following this procedure, 

however, one must be aware of the following subtlety in the equations. 

For each q = j + i the set of equations 



( |B9| ) with 1 = 1 and j = q — 1, and 
( Bill) for all i = 1,3, . . . ,q with j = q — i 



can be shown to be linearly dependent for arbitrary k and for any equilibrium stellar model. For 
example, taking the simplest case of g = 1, one can show that equation ( |B9| ) with i = 1 and j = 
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becomes 



= (m + 2) [wi^i - (m + l)vi^i] 



while equation (Bll) with i = 1 and j = becomes 



= Qm+iCm_2|^K(m + 2)cmCm+2 +m[{m + 1)0^ + bm] - (2m + 3)(m + 3)'^Ql^^2Cm I 
X [wi^i - (m + . 

This problem can be solved by eliminating one of these equations from the subset for each q 



(for example, equation (|BllD with i = 1). Thus, to properly construct the algebraic system ( [46| ) 
we use, for all j = 0, 2, . . . , j^^^, the equations 



(EH) 
(H) 



( |B9| ) with i = 1, 3, 

(plOl) with A; = 0,1 

(pill) with i = 3,5, 

(pl2|) with A; = 0,1 



; * * * ! '"max 



; * * * ! '"max 



B.2. Polar Hybrids 

For I = m, m + 2, m + 4, . . . the regular series expansions'^ about the center of the star, r = 0, are 



Wm+jir) 
Um+j+i{r) 



R 



m+j 



r) 



1=0 
i even 

OO 

E 

i=0 
i even 

OO 



i even 



Rj 



(B13a) 
(B13b) 
(B13c) 



where j = 0, 2, 4, . . .. 

The regular series expansions about r = R, which satisfy the boundary condition Ap = are 



Wm+j{r) 



|%'.(i-^)' 



(B14a) 



^^We present the form of the series expansions for Ui{r) for reference; however, we do not need these series since 
we eliminate the Ui{r) using equation (pSl). 
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k=0 

oo 



k=0 



T \ ^ 



R 



(B14b) 
(B14c) 



where j = 0, 2, 4, . . .. 



These series expansions must agree in the interior of the star. We impose the matching condition 
that the series ( [B13 ) truncated at i,„ax be equal at the point r = tq to the corresponding series 



( B14 ) truncated at k^^^^. That is, 



- J2 - 



R 



z even 



R J ^ ' ^ \ R / 7 n 

i even 



R 



(B15a) 
(B15b) 



When we substitute (^T|) and (||) into (H), the coefficient of (r/i?)'"+^+* in the resulting 
equation is 



i-2 



= (m + j + i + l)wj^i + ^ TTs+i Wj^i-s-2 - (m + j)(m + j + l)i 



(B16) 



Similarly, when we substitute ( B14| ) and ( p5[ ) into (^), the coefficient of [1 — {r/R)]^ in the 
resulting equation is 

k 

= (A; + 1) [wj^k - Wj,k+i\ + X! i^s-i - Wj^k-s+i - {m + j){m + j + l)vj^k (B17) 

s=0 

where we have defined -7f_2 = = wjfi. 

When we use (^) to eliminate the Ui{r) from ( |39| ) and then substitute for the Wi±i{r) and Vi±i{r) 
using ( [B13| ), the coefficient of (r/i?)™'^-'+* in the resulting equation is 



-im{m + j - l)Qm+jQm+j~lCm+j+l 



Wj-2,i+2 - {m + j - 2)vj-2,i+2 



+ \{m + j - l)Ql^^j[{i + l)m - iK(m + j - l){m + j + i)]cm+j+i 



Cm+j — 1 Cm+j + 1 



7i + j + i){l- Q^+j- - Qm+j+i) + ^'^f 
- (m + i + 2)Q^+j.+i[m(2m + 2j + i + 2) + \K{m + j + 2)(m + j + i)]c„+j_i|u;j,i 



+ |(m + j - l)(m + j + l)Ql,^j[{i + l)m - ^K{m + j - l){m + j + i)]cm+j+i 

- rn^ + {m + j + i)am+j Cm+j-iCm+j+i 
+ (m + i)(m + j + 2)Q^+^.+i 

X m{2m + 2j + i + 2) + + j + 2)(m + j + z 

+ (5r„+j+2(3m+i+i + j + i)+ m{m + j + l)(2m + 2 j + i + 2)] Cm+j-i 

Wj+2,i-2 + (m + j + 3)T;j+2,i-2 



(B18) 



When we use (^) to eliminate the Ui{r) from ( p9| ) and then substitute for the Wi±i{r) and Vi±i{r) 
using (|B14|), the coefficient of [1 — (r/R)]'' in the resulting equation is 



= m{m + j - l){m+ j - k)Qm+jQm+j~iCm+j+i 



+ {k + l)m{m+j - l)Qm+jQm+j-lCm+j+l 



Wj-2,k+i - {m + j - 2)vj_2, 



fc+i 



+ < -{m + j - l)Q'^^j[{^Kk + m){m + j - 1) - km]cm+j+i 



+ 



'^nm + 



^ Qm+j Q 



Cm+j — 1 ^m+j + 1 



{m + j + 2)Q^+j+i [{^Kk + m){m + j + 2) + k 



m\Cm+j-l (Wj^k 



- (A; + + j - l)Ql^^j[m - lK{m + j - l)]cm+j+i 

+ ^1 — Qm+j ~ Qm+j+lj Cm+i-lCm+jr+l 

- {m + j + 2)Ql,_^j^^[m + lK{m + j + 2)]cm+j-i^Wj^k+i 

+ |-(m + j - l){m + j + l)Q'^^j[{^Kk + m){m + j - 1) - km]cm+j+i 
— {m? + ka^+j^ Cm+j-iCm+j+i 

+ {m + j){m+j + 2)Q^+j+i[(iKA: + m)(m + j + 2) + km]cm+j-i^Vj^k 
+ {k + l)\^-{m + j - l){m + j + l)Q^+j[m - ^K{m + j - l)]cm+j+i 

+ Oim~{-jCm+j~l(^m+j+l 
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(m + j){m +j + 2)Q.^+j+i [m + ^/{(m + j + 2)]cm+j_i jiij- fc+i 



+ m{m + j + 2)(m + j + k + l)Qm+j+2Qm+j+ic, 



m+i — 1 



{k + l)m{m + j + 2)Qm+j+2Qm+j+ic, 



m+j — 1 



W 



'j+2,k+i + {m + j + 3)vj+2,k+i 



(B19) 



The equations ( B15| ) through ( |B19| ) make up the algebraic system ( p6| ) for eigenvalues of the 
polar-led hybrid modes. As in the case of the axial-led hybrids, one truncates the angular and 
radial series expansions at indices jmaxj ^max and /c^ax and constructs the matrix A by keeping the 
appropriate number of equations for the number of unknown coefficients Wj^i, vj^i, Wj^k and Vj^k- 

We, again, find that certain subsets of these equations are linearly dependent for arbitrary n 



and for any equilibrium stellar model. For all j, it can be shown that both equation ( B16 ) with 



i = and equation ( B18 ) with i = are proportional to 



This problem can, again, be solved by eliminating, for example, equation ( [B18| ) with i = for 
all j. Thus, to properly construct the algebraic system (|^) we use, for all j = 0,2, . . . , jmax, the 
equations 

dmi) 

(|BT5b| ) 

(IbTgI) with i = 0, 2, 



-';***! "max 



(|BT7|) with A; = 0, 1, , 



(|BT8| ) with i = 2, 4 



? • • • ) "max 



(|BT9| ) with A: = 0,1,...,A;^,. -1. 
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Fig. 1. — All of the non-zero coefficients W/(r), V/(r), Ui{r) of the spherical harmonic expansion 
( p^ ) for a particular m = 2 axial-led hybrid mode of the uniform density star. The mode has 
eigenvalue kq = —0.763337. Note that the largest value of / that appears in the expansion ( p3| ) 
is = 4 and that the functions Wi{r), Vi{r) and C//(r) are low order polynomials in (r/R). (See 
Table ^.) The mode is normalized so that V2{r = R) = 1. 

Fig. 2. — The coefficients Wi{r), V/(r), Ui{r) with / < 6 of the spherical harmonic expansion ( p3|) 
for a particular m = 2 axial-led hybrid mode of the polytropic star. This is the polytrope mode 
that corresponds to the uniform density mode displayed in Figure For the polytrope the mode 
has eigenvalue kq = —1.025883. The expansion (^) converges rapidly with increasing I and is 
dominated by the terms with 2 < ^ < 4, i.e., by the terms corresponding to those which are non- 
zero for the uniform density mode. Observe that the coefficients shown with 4 < / < 6 are an order 
of magnitude smaller than those with 2 < / < 4. Those with I > 6 are smaller still and are not 
displayed here. The mode is, again, normalized so that V2{r = R) = 1. 

Fig. 3. — The functions Wi{r) with / < 6 for a particular m = 1 polar-led hybrid mode. For 
the uniform density star this mode has eigenvalue kq = 1.509941 and Wi = —x + (x = r/R) 
is the only non- vanishing Wi{r) (see Table §). The corresponding mode of the polytropic star 
has eigenvalue kq = 1.412999. Observe that Wi{r) for the polytrope, which has been constructed 
from its power series expansions about r = and r = R, is similar, though not identical, to the 
corresponding Wi{r) for the uniform density star. Observe also that the functions Wi{r) with / > 1 
for the polytrope are more than an order of magnitude smaller than Wi{r) and become smaller 
with increasing I. (VF5(r) is virtually indistinguishable from the {r/R) axis.) 

Fig. 4. — The functions Vi{r) with / < 6 for the same mode as in Figure ^. 

Fig. 5. — The functions Ui{r) with / < 6 for the same mode as in Figure ^. 

Fig. 6. — The functions Ui{r) with I < 7 for a particular m = 2 axial-led hybrid mode. For the 
uniform density star this mode has eigenvalue kq = 0.466901 and U2{r) and U4,{r) are the only non- 
vanishing Ui{r). (See Table ^ for their explicit forms.) The corresponding mode of the polytropic 
star has eigenvalue kq = 0.517337. Observe that U2{r) and Uii{r) for the polytrope, which have 
been constructed from their power series expansions about r = and r = R, are similar, though 
not identical, to the corresponding functions for the uniform density star. Observe also that the 
UQ{r) is more than an order of magnitude smaller than U2{r) and U/^{r). 

Fig. 7. — The functions Wi{r) with / < 7 for the same mode as in Figure |^. 

Fig. 8. — The functions Vi{r) with / < 7 for the same mode as in Figure ^. 

Fig. 9. — The functions Ui{r) with / < 8 for a particular m = 2 axial-led hybrid mode. For 
the uniform density star this mode has eigenvalue kq = 0.359536 and U2{r), U4{r) and UQ^r) are 
the only non- vanishing Ui{r). (See Table ^ for their explicit forms.) The corresponding mode 
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of the polytropic star has eigenvalue kq = 0.421678. Observe that U2{r), U4{r) and (/^{r) for 
the polytrope, which have been constructed from their power series expansions about r = and 
r = R, are similar, though not identical, to the corresponding functions for the uniform density 
star. Observe also that Us{r) is more than an order of magnitude smaller than U2{r), U4{r) and 
Ueir). 

Fig. 10. — The functions Wi{r) with / < 8 for the same mode as in Figure ^. 
Fig. 11. — The functions Vi{r) with / < 8 for the same mode as in Figure 
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6253:^ 




11.253:4 _|_ g 


6253:6 





623029(53:4 


- 73:6) 




1.861684 





778500(8.753:^ — 31.53:"^ + 24. 75^^) 


-0.326313(9:^4 _ ^^^6-^ 





168645^^ 


5 


25^2 


— 13.53:4 _|_ g 


.253:^ 


-0 


192134(203:4 


- 203:6) 


2 


6253:^ 




11.25^4 _|. 9 


6253:6 


-0 


192134(5^4 


- 73:6) 


^The 


TTi + 1 = io i 


9 the 


normalized so that V2 = 1 at the su 
maximum value of i in the spherica 


rfaee of the star, ai = 1 . He 
harmonic expansion of 


Obse 


T./R). 

ve that this 


i = 


^ Iq term is always 



























TABLE 2 

Polar-Hybrid Eigenfunctions'^ with m = 1 for Uniform Density Stars. 





^■^0 








1 












V3(r) 








U4(r) 


2 -0, 


,176607 






-a: + 3 


j,3 











-X + 2a;3 







-0.876991.^ 







1, 


,509941 






+ 3 


„3 











-X + 2a;3 







0.380087.^ 







4 -1, 


,183406 


1, 


,25a3 


- 3.5a3^ 


+ 2.25a3^ 


-1.585327(6^^ 


- 6a;^) 


1, 


,25a3 


- + 6.75a;^ 


-1.585327(2a;3 


- 3.^) 


0.813707(7.^ - 9.^) 


-0 


.90411033^ 


-0. 


.068189 


1. 


.25a; 


- 3.5a;3 


+ 2.26a;^ 


0.100030(6a;3 


- 6a;5) 


1. 


.25a; 


- 7a;3 + 6.75a;^ 


0.100030(2a;3 


- 3.^) 


0.398091(7.3 - 9.^) 





.435309.^ 


1. 


.045597 


1. 


.25a; 


- 3.5a^3 


+ 2.25a^'^ 


0.331793(6.3 


- 6a;5) 


1. 


.25a; 


- 7a;3 + 6.75a;S 


0.331793(2.3 


- 3.''^) 


-0.010993(7.3 - 9.''^) 


-0 


.250819.^ 


1, 


,805998 


1, 


,25a; 


- 3.5.3 


+ 2.25a;^ 


-0.343160(6a;3 


- 6a;5) 


1, 


,25a; 


- 7a;3 + 6.75a;^ 


-0.343160(2.3 


- 3.'-^) 


-0.300378(7.3 _ q^5j 





.147226.^ 



The eigenfunctions are normalized so that = 1 at the surface of the star, . = 1. Here . = (r/R). 

^Iq — m -\- 1 — Iq is the maximum value of I in the spherical harmonic expansion of Sv^ . Observe that this I — Iq term is always axial. 



TABLE 3 

AxTAT -H^TiTiin ET<"irxFT"x<"'TTOxs" AVTTTT ))) — 2 r<")Ti TJxTr<")Tn[ DrxsTT^- Stars. 



00 



do - 1)^ 




«0 










U4(-r) 






Uq (-r) 






W3(r 








W5(r-) 






V3(r) 








1 





666667 




^3 














































3 


-0 


763337 




0.352414(73^^ - Om^) 






-0.679569£c^ 













6(.3 - 


3:5) 













— 23:3 _|. 3^,5 














466901* 




2.522714(7iP^ — 0=^^ ) 






2.452800^:^ 













6(^.^3 _ 















— 2ip3 _|_ 3^.5 











1 


496436 




-0.607406(7£c^ - Om^) 






0.5049643:^ 













6(.3 - 















— 23:3 _|_ 3^5 









5 


-1 


308000 


-0 


510418(7. 875£c2 - 24.75=:^ + 17 


8753:''') 





634277(113:^ - 


13.7) 


-0 


6396093:7 


13 


5.3 


- 33.^ 


+ 19 


5.7 


- 1 


138387(15.5 


- 15.7) 




53:3 _ ie.5a:5 + 133,7 




138387(33:5 


- 4.7) 




-0 


GQggg4 


-0 


SB6BSl(7.S75a:^ — 2-! . Tr>.r + 17 


875^''') 





188642(lla:^ — 


13.7, 





4558273:7 


13 


5.3 




+ 19 


5.7 





349918(16.5 


- 15.7) 




r-..-'--'- — ir-,.r-.s:° + 133:7 





349918(33:5 


- 4.7) 







3B9B36* 


-3 


28 1 6 79(7.8 75a: — 2-1.T.-..r'^ + 17 


875^''') 


-0 


769879(lla:^ — 


13.7) 


- 1 


1038003:7 


13 


5.3 




-) 19 


5.7 





37Q022(1B.5 


- 15.7) 




- ir>.,-,a:5 + 133:7 





370022(33:5 


- 4.7) 




1 


1B30BS* 


2 


072212(7. 875ic2 — 24.75a:^ + 17 


875^''') 





102912(113!^ — 


13.7) 


-0 


5736893:7 


13 


5.3 




+ 10 


5.7 





769719(lB!n5 


- 15.7) 




5a:3 _ 16.53:5 ^ I3a.7 





769719(3a:5 


- 4.7) 




1 


733971* 





944346(7. 875ic2 — 24.75ic^ + 17 


875^''') 


-0 


423766(ll!B^ — 


13.7) 





280929a:7 


13 


5.3 


- 33.^ 


+ 19 


5.7 


-0 


583914(lB!n5 


- 15.7) 




53:3 _ ig.5a,5 ^ 13a,7 


— 


583914(3!b5 


- 4.7^ 



°-The eigenfunctions are normalized so that V3 = 1 at the surface of the star, a: = 1 . Here at = (r / R) . 

^ Iq — m + 1 ^ Iq — 1, where Iq is the maximum value of I in the spherical harmonic expansion of Su" . Observe that this i 



TABLE 4 

Polar-Hybrid Eigenfunctions'^ with m = 2 for Uniform Density Stars. 







KQ 






W2(r) 




W4(r) 








V.2(r) 






V4(r) 




U3(r) 






2 


-0 


.231925 






-3.2 _^ 3.4 













1.5.2 2. 


5.4 









-0.891544.4 









1 


.231925* 






-3.2+3.4 













1.5.2 _|_ 2_ 


5.4 









0.560825.4 







4 


-1 


.092568 


5, 




- 13.5.4 _|_ 5 


i.25.^ 


-0.909581(20.4 


- 20.6) 


2, 


,625.2 


- 11.25.4 


+ 9 


.625.6 


-0.909581(5.4 


- 7.6) 


0.872718(9.4 _ ii^6-, 


-1 


.093523.6 




-0 


.101790 


6. 




- 13.6.4 + S 


i.25.^ 


0.078913(20.4 


- 20.6) 


2. 


.625.2 


- 11.25.4 


+ 9 


.625.6 


0.078913(6.4 


- 7.6) 


0.381215(9.4 - 11.6) 





.494643.6 




0. 


.884249* 


6. 


.25x^ 


- 13.6.4 + S 


i.25.^ 


0.176440(20.4 


- 20. <3) 


2. 


.625.2 


- 11.25.4 


+ 9 


.625.'^ 


0.176440(6.4 


- 7.6) 


-0.107938(9.4 - 11.6) 


-0 


.346296.6 




1 


.643443* 


5, 


.25x^ 


- 13.5.4 _|_ 5 


i.25.^ 


-0.350886(20.4 


- 20. <^) 


2, 


,625.2 


- 11.25.4 


+ 9 


.625.6 


-0.350886(5.4 


- 7.6) 


-0.484558(9.4 _ 113,6) 





.342451.6 



The eigenfunctions are normalized so that V2 = 1 at the surface of the star, . = 1. Here . = (r / Ft) . 
^Iq — m -\- 1 — Iq — 1, where Iq is the maximum value of I in the spherical harmonic expansion of Sv^ . Observe that this I — Iq term is always axial. 
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Table 5. Eigenvalues kq^ for Uniform Density Stars. 



(lo - m + 1)'' 


parity'^ 


m = 


m = 1 


m = 2 


m = 3 


m = 4 


1^ 


a 


0.000000 


1.000000 


0.666667* 


0.500000* 


0.400000* 


2 


P 


-0.894427 


-0.176607 


-0.231925 


-0.253197 


-0.261255 




P 


0.894427 


1.509941 


1.231925* 


1.053197* 


0.927922* 


3 


a 


-1.309307 


-0.820009 


-0.763337 


-0.718066 


-0.680693 




a 


0.000000 


0.611985* 


0.466901* 


0.377861* 


0.317496* 




a 


1.309307 


1.708024 


1.496436* 


1.340205* 


1.220340* 


4 


P 


-1.530111 


-1.183406 


-1.092568 


-1.022179 


-0.965177 




P 


-0.570463 


-0.068189 


-0.101790 


-0.120347 


-0.131215 




P 


0.570463 


1.045597 


0.884249* 


0.773460* 


0.691976* 




P 


1.530111 


1.805998 


1.643443* 


1.511923* 


1.404416* 


5 


a 


-1.660448 


-1.404217 


-1.308000 


-1.230884 


-1.167037 




a 


-0.937698 


-0.537334 


-0.509994 


-0.486868 


-0.466934 




a 


0.000000 


0.440454* 


0.359536* 


0.304044* 


0.263530* 




a 


0.937698 


1.306079 


1.153058* 


1.040073* 


0.952507* 




a 


1.660448 


1.861684 


1.733971* 


1.623634* 


1.529045* 



^kqQ, — {a + mO.) is the mode frequency in the rotating frame to lowest order in Q. The 
modes whose frequencies are marked with a * satisfy the condition a{a + mO.) < and are 
subject to a gravitational radiation driven instability in the absence of viscous dissipation. 

^For m = 0, this is simply Iq. For the uniform density star, lo is the maximum value of I 
appearing in the spherical harmonic expansion of Sv"^ . 

'^This denotes the parity class of the mode; a meaning axial-led hybrids, and p meaning 
polar- led hybrids. 

'^These are the eigenvalues of the pure lo = m r-modes. For isentropic stars they are 
independent of the equation of state and have the value kq = 2/(m -I- 1) (or kq = for 
m = 0) to lowest order in SI ( Papalouizou and Pringle 197^ ). 
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Table 6. Eigenvalues kq'^ for the p = Kp^ Polytrope. 





parity'^ 


m = 


m = 1 


m = 2 


m = 3 


m = 4 


1" 


a 


0.000000 


1.000000 


0.666667* 


0.500000* 


0.400000* 


2 


P 


-1.028189 


-0.401371 


-0.556592 


-0.631637 


-0.672385 




P 


1.028189 


1.412999 


1.100026* 


0.904910* 


0.771078* 


3 


a 


-1.358128 


-1.032380 


-1.025883 


-1.014866 


-1.002175 




a 


0.000000 


0.690586* 


0.517337* 


0.412646* 


0.342817* 




a 


1.358128 


1.613725 


1.357781* 


1.176745* 


1.041683* 


4 


P 


-1.542065 


-1.312267 


-1.272885 


-1.238631 


-1.208390 




P 


-0.701821 


-0.178792 


-0.275335 


-0.333267 


-0.370450 




P 


0.701821 


1.051525 


0.862948* 


0.734297* 


0.640592* 




P 


1.542065 


1.726257 


1.519573* 


1.360560* 


1.234698* 


5 


a 


-1.656481 


-1.483402 


-1.433916 


-1.391305 


-1.354057 




a 


-1.013703 


-0.705182 


-0.703898 


-0.699942 


-0.694498 




a 


0.000000 


0.528102* 


0.421678* 


0.350192* 


0.299055* 




a 


1.013703 


1.281962 


1.104402* 


0.974192* 


0.874124* 




a 


1.656481 


1.795734 


1.627215* 


1.489441* 


1.375406* 



'^kqQ — {a + mO.) is the mode frequency in the rotating frame to lowest order in Q. The 
modes whose frequencies are marked with a * satisfy the condition a{a + mil) < and are 
subject to a gravitational radiation driven instability in the absence of viscous dissipation. 

''For m = 0, this is simply Iq. For the n = 1 polytrope, lo is the largest value of I that 
contributes a dominant term to the spherical harmonic expansion of 5v°' . 

'^This denotes the parity class of the mode; a meaning axial-led hybrids, and p meaning 
polar- led hybrids. 

'^These are the eigenvalues of the pure lo = m r-modes. For isentropic stars they are 
independent of the equation of state and have the value kq — 2/(m -I- 1) (or kq = for 
m = 0) to lowest order in Q (Papalouizou and Pringie 1978). 
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Table 7. Dissipative timescales (in seconds) for m = 1 axial-hybrid modes^ at T = l(fK and 



lo 








fs 






3 





0.611985 




7.67 X 10^ 


-9.79 X 10® 






1 


0.690586 


5.86 X 10^ 


9.29 X lO'' 


-1.25 X 10® 


-1.22 X 10^° 


5 





0.440454 




2.04 X 10^ 


— oo 


-2.07 X 10^^ 




1 


0.528102 


2.57 X 10^ 


3.87 X lO'' 


-2.17 X 10^° 


-5.75 X 10" 



^We present dissipative timescales only for those modes that arc unstable to 
gravitational radiation reaction. None of the m = 1 polar-hybrid modes are 
unstable for low values of Iq. 

''The polytropic index, n, where p = Kp^'^^^". The n=0 case represents the 
uniform density equilibrium star. 

'^Dissipation due to bulk viscosity is not meaningful for uniform density stars. 
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Table 8. Dissipative timescales (in seconds) for m = 2 axial-hybrid modes at T = lO'^A' and 



{lo - 1) 


n" 


K 


Tb 


fs 




f2 


f4 


re 









0.666667 




4.46 


X 


10* 


-1.56 X 10° 












1 


0.666667 


2.0 X 10" 


2.52 


X 


10® 


-3.26 X 10° 










3 





0.466901 




4.10 


X 


10^ 


— oo 


-3.88 X 10^ 










1 


0.517337 


6.43 X 10^ 


6.21 


X 


10^ 


< -10^® 


-1.85 X 10® 


-4.97 


X 


lO^s 







1.496436 




3.92 


X 


10^ 


—oo 


-5.85 X 10^ 










1 


1.357781 


4.10 X 10^ 


7.18 


X 


10^ 


< -10i« 


-1.60 X 10^ 


-4.35 


X 


10i« 


5 





0.359536 




1.34 


X 


10^ 


— oo 


— oo 


-1.28 


X 


10" 




1 


0.421678 


2.65 X 10^ 


3.01 


X 


10^ 


< -lO^^ 


-2.01 X 10^ 


-1.15 


X 


10^2 







1.153058 




1.32 


X 


10^ 


— oo 


— oo 


-3.11 


X 


10" 




1 


1.104402 


2.45 X 10^ 


3.65 


X 


10^ 


< -10^^ 


-1.37 X 10" 


-4.89 


X 


10" 







1.733971 




1.31 


X 


10^ 


— oo 


— oo 


-1.92 


X 


10^1 




1 


1.627215 


5.32 X 10^ 


3.44 


X 


10^ 


< -10i« 


-2.30 X 10^^ 


-8.33 


X 


10i« 



^The polytropic index, n, where p = Kp^^^^^. The n=0 case represents the uniform density 
equilibrium star. 

''Dissipation due to bulk viscosity is not meaningful for uniform density stars. 

"^This is the lo — m = 2 r-mode aheady studied by Lindblom et al. (1998), Owen ct al. (1998), 
Andersson et al. (1998), Kokkotas and Stergioulas (1998) and Lindblom et al. (1999). The value 
of the bulk viscosity timescale for this mode is taken from Lindblom et al. (1999) who calculate it 
self-consistently using an order Q,^ calculation. 
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Table 9. Dissipative timescales (in seconds) for m = 2 polar-hybrid modes at T = lO^K and 



(Zo - 1) 






fa" 








T5 


2 





1.231925 




9.03 X 


10^ 


-4.77 X 10'' 






1 


1.100026 


3.32 X 10^ 


1.24 X 


10* 


-3.37 X 10'' 


-3.13 X 10" 


4 





0.884249 




2.17 X 


10^ 


— oo 


-5.64 X 10^ 




1 


0.862948 


1.93 X 10^ 


4.94 X 


10^ 


-1.10 X 10^ 


-1.45 X 10^° 







1.643443 




2.13 X 


10^ 


— oo 


-2.12 X 10^^ 




1 


1.519573 


4.79 X 10^ 


4.77 X 


10^ 


-1.92 X 10" 


-2.31 X 10" 



'^The polytropic index, n, where p = Kp^~^^^". The n=0 case represents the uniform 
density equilibrium star. 

''Dissipation due to bulk viscosity is not meaningful for uniform density stars. 




Figure 1 




Figure 2 



-48- 



0.5 



G 


O 







Wi{r) for the uniform density star 
Wi (r) for the n = 1 polytrope 
W3{r) for the n = 1 polytrope 
Ws{r) for the n = 1 polytrope 



C^- 



-0.5 







0.5 

(r/R) 



Figure 3 



-49- 




-0.5 ' ' ' ' ' ' ' ' ' ' ' 

0.2 0.4 0.6 0.8 1 

(r/R) 



Figure 4 
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Figure 6 
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Figure 7 
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Figure 8 
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Figure 9 
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Figure 10 
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